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Abstract 

We present the results of the U BV RqIq surface brightness profiles de- 
composition of the Seyfert galaxies III Zw 2, Mrk 506 and Mrk 509. The 
profiles were modelled as a sum of a Gaussian, a Sersic law and an expo- 
nent. A Ferrers bar and a Gaussian ring were added to the model profiles of 
III Zw 2 and Mrk 506, respectively. The parameters and the total magnitudes 
of the structural components were derived. 



1 Introduction 

We present here the first results of an outgoing study aimed to do a detailed de- 
composition of Seyfert galaxies UBVRcJc surface brightness profiles (hereafter 
SBPs). Our advantages over Seyfert galaxies SBP decomposition found in the 
literature are the following: (1) we model explicitly the active nucleus in oppo- 
sition to some authors who avoid nucleus modelling (e.g. Q), (2) we use Sersic 

*Based on observations obtained at the Rozhen National Astronomical Observatory of Bulgaria 
operated by the Institute of Astronomy, Bulgarian Academy of Sciences. 
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rather than de Vaucouleurs law for bulge, (3) we use a truncated exponential law 
defined in lEl along with a pure exponent in disk modelling, and (4) we model 
bar/oval/lens/ring components that have been generally skipped by Seyfert galax- 
ies SBPs decomposers (e.g. [3|). 

All galaxies to be decomposed were observed at Rozhen NAO of Bulgaria 
with the 2-m telescope and Photometries AT200 CCD camera (0.309 arcsec px"^) 
through a standard Johnson-Cousins UBVRcIq set of filters. The SBPs of the 
galaxies were extracted fitting ellipses to the galaxian isphotes by means of fit/ells 
command of the surfphot context of eso-midas package (see [4| for details). 

Our first decomposition results concern the Seyfert galaxies III Zw 2 (Syl.O), 
Mrk 506 (Syl.5) and Mrk 509 (Syl.2); the Seyfert types are taken from NED. 

2 Models and Methods 

Our basic model SBP is a sum of (1) a Gaussian with a fixed FWHM to repre- 
sent the nucleus (the only free parameter is the central surface brightness), (2) a 
Sersic law [51 with free parameters /ics - the effective surface brightness, r^f^ - 
the effective radius, and n - the power-law index, to represent the bulge and (3) 
an exponent ^ with free parameters /icen - the central surface brightness, and 
Tsci - the scale length, to represent the disk. Bulge and disk model SBPs were 
convolved with a Gaussian PSF to simulate the seeing effects on the profiles ac- 
cording to fT'JHl. Note that if the frame PSF is not circular then the FWHM of the 
convolution Gaussian is set to the PSF's mean FWHM along the minor axis and 
the FWHM of the nuclear Gaussian is set to the PSF's mean FWHM along the 
major axis (the mean FWHM of the stellar images for each frame was determined 
fitting a 2D Gaussian to a number of field stars employing center/iqe command 
within eso-midas). 

We added (1) a Ferrers model profile [9 1 with free parameters /icen - the central 
surface brightness, rend - the profile length, and m - the power-law index, to 
account for the bar^ in III Zw 2, and (2) a displaced Gaussian model profile [T0| 
with free parameters /icen - the central surface brightness, rcen - the position (or 
displacement) of the Gaussian centre, and JTW - the FWHM, to account for the 
ring in Mrk 506 (the galaxy has SAB(r)a morphology according to NED). 

The figure-of-merit function that is minimized is equal to the unweighted sum 
of the squared differences between the observed and the model SBPs per degree 

'The bar manifests in about 0.25 rise of the elHpticity and an almost constant position angle. 
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of freedom, v: 

N 



i=l 



(1) 



where p is the P-element vector of the free model parameters to be fitted. The 
initial value of the degrees of freedom is defined d& uq = N — P, where N and 
P are the number of the profile data points and the number of the free parame- 
ters, /?, to be fitted, respectively. Note that corrections to vq could be done for 
the presence of (1) zero-weighted profile data points, A^zcro, and/or (2) fixed pa- 
rameters, Pfixj so the actual value of v becomes z/ = z/q — A^zoro + -Pfix. and the 
corrected value of v enters Eq.[T] The minimization of A^(/7) was performed em- 
ploying Levenberg-Marquardt algorithm fTl\. The initial guess parameters were 
estimated by eye overplotting the observed and the model profiles and changing 
the parameters manually to get a good correspondence between them. After that, 
a P-dimensional valley \p — 6p ,p + 6p~^] around the initial guess parameters was 
constructed; here 5p^/^ are the allowed parameters' deviations from their initial 
guess values in both directions (these deviations could be different for each param- 
eter). Next, a number of decomposition cycles were run: in each cycle the actual 
initial guess parameters were picked up randomly from the uniformly distributed 
parameters in the valley defined above. This procedure helped us in isolating the 
global minimum among /^lip) values - the larger the number of the random cy- 
cles is, the bigger the probability to find the global minimum becomes. The num- 
ber of random cycles could vary depending on the complexity of the profile to 
be decomposed - more complicated profiles could require up to several hundreds 
of random cycles. After the minimization cycles were completed a histogram of 
^lip) values was built and the minimum corresponding to the most frequently 
occurring /^lip) was selected. If the parameters corresponding to this minimum 
behaved themselves well and if this minimum had the lowest value of A^(j9) then 
it was assumed to be the global one with Al(p^^^^) = the correspondence 

between the parameter values obtained after decomposition of different observing 
runs profiles could be used as a further check of the minimum found (see Table[Tl). 
If some of the parameters corresponding to the most frequently occurring mini- 
mum had unacceptable values (e.g. very small or very large), then this minimum 
was rejected and a new global minimum was searched for among the remaining 
random cycles results. If there were minima with A^(/>) values lower than A^(p) 
of the most frequently occurring minimum then these minima were checked one 
by one; in all cases we found that the parameters corresponding to these minima 
had unacceptable values, i.e. these minima were local ones. 
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3 Results 



At the end of the selected objects decomposition the best-fit parameters were ob- 
tained (listed in Tabled]) and, based on them, the total magnitudes of the structural 
components were computed (listed in Table EJ. We list only the total magnitudes 
for the nuclear Gaussian because its central surface brightness is strongly depen- 
dent on the seeing and does not allow straightforward comparison of the results 
obtained at nights with different seeing conditions. The structural parameter(s) 
and/or the total magnitude(s) that could not be derived from the decomposition 
are marked in the tables. The errors of the parameters are la uncertainties as re- 
sulting from the fitting algorithm and they should be considered as approximate 
ones. More reliable estimate of the parameter errors could be obtained through 
Monte Carlo simulations or a bootstrap analysis. The mean FWHM along the 
minor axis of the stellar images, JFVVpsf, that was used in the model SBPs con- 
volution and the values of (jfit = (A^ min)°'^ listed in Table |2l as well. The 
civil date and the name of the decomposed object are listed in the tables using 
the following code: la - September 9/10, 1997 III Zw 2; lb - June 1/2, 1997 
Mrk 506; 2b - July 18/19, 1998 Mrk 506; Ic - July 10/11, 1997 Mrk 509; 2c - 
September 8/9, 1997 Mrk 509; 3c - July 20/21, 1998 Mrk 509. Corresponding 
Johnson-Cousins filter is shown along with the date and the object code. We list 
the parameters of the bar in III Zw 2 and of the ring in Mrk 506 in Table |3] and 
TablelH respectively. Note that the magnitudes and the surface brightnesses listed 
in the tables have not been corrected for the Galactic absorption and cosmological 
dimming; k- and evolution corrections have not been applied as well. 

The observed SBPs and the decomposed profiles of the structural components 
of m Zw 2, Mrk 506 and Mrk 509 are shown in Fig. □ Fig. El and Fig. re- 
spectively, along with the residual profiles equal to — fJ'^°°ip rnin)- We show 
the decompositions with the smallest a^n among the different filters and observing 
runs for each galaxy (see Table |2l). 

We list in Table|3]the length of the bar in III Zw 2 obtained from the EVRqIc 
ellipticity profiles using ( 1 ) a maximum ellipticity criterion - the bar length, /^^'^''^ , 
corresponds to the point of maximal ellipticity, and (2) a minimum ellipticity cri- 
terion - the bar length, /b^""*, corresponds to the point of minimal ellipticity next 
to the ellipticity maximum, so one has l^^^^^^ > ^b™'^''^ definition. One could 
see that /b™™'* agree well with the bar length, rend, obtained from the SBPs de- 
composition, while /b™^'''' underestimates rend (cf- [T^)- 

BVIc profiles decomposition of Mrk 509 was recently presented by where 
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Table 1: Bulge (°) and disk (°) structural parameters for III Zw 2, Mrk 506 and 
Mrk 509. The errors of the parameters are shown as a superscript index to the 
corresponding values. 

Code ^i^s r^s n« r^^, 

[mag/D"] ["] [mag/D"] ["] 

M 18.7620-193 0.8490058 0.976 o-^^^ 23.3800-028 5.i380-044 

UV 18.1820-143 0.8470-047 1.1170.104 22.00l0-oi7 5.732O-033 

lai?c 17.3460-315 o.7840-092 0.915O-I86 21.1460-05i 5.3180-085 

la/c 17.9390-340 1.1140-138 0.6I5O-177 20.3120023 5.766O-042 



2hU 


18.1810-210 


1.0030-075 


1.0790-140 


21.4040053 


6.4730-124 


IhB 








21.0110-010 


6.2220-021 


2hB 


19.8250-242 


1.3420-153 


0.3160-171 


21.1100-023 


6.1850-045 


IhV 








20.2270-007 


6.3460-013 


2hV 


19.1680-107 


1.2920-054 


0.5920-077 


20.1720.016 


6.1980-031 


IhRc 


19.563^-958 


1 2871-324 


1.1921-694 


19.5790013 


6.0510-022 


2hRc 


19.0400-043 


1.473 0-029 


0.4600-036 


19.6280-010 


6.2850-018 


Ib/c 


19.6 1 90-050 


2.0570-054 


0.6470-083 


19.1610-009 


7.0860-023 


2b7c 


18. 149 0-034 


1.2980-019 


0.6580-036 


18.9390-014 


6.2710-030 



3cU I6.9150-201 1.5450-062 1.2030-042 

IcB 18.156 0-086 1.4290-045 1.3590-065 21.9330-045 5.7630-092 

2cB 17.1840-091 0.9810-032 i.8oo0-063 21.84800i5 6.4250-033 

3cB 17.4170-212 1.45 0-074 1.108 0-044 2I.4IIO-022 5.546O.042 

IcV I8.4430-066 1.88 0-047 0.85 0-044 2O.514O-018 5.046O-025 

2cV 18.2 1 00-082 1.58 5 0-049 i.3300-075 20.9100-027 6.6540-070 

3cV 17.4140-144 1.4000-070 i.5080096 20.8770-028 6.253O-O60 

lci?c 17.7650-034 1.7170.023 1.0530-028 20.0510011 5.3060-oi7 

2ci?c 17.7440-049 1.6330-032 1.4050-049 20.5500-oi6 6.841 oo35 

3ci?c 17.88 6 0-072 1.8 7 7 0-048 i.02l 0-048 20.10 7 0-oi8 5.9730-039 

Ic/c 17.63 2 0-030 1.9570.024 0.84 1 0-028 19.217O-013 5.245 0-022 

2c/c 17.9770-045 1.96 0-034 1.08 3 0-042 19.6I4O-011 6.4430-024 

3c7c 17.7510-057 2.0360-044 0.951 o-049 19.356O.018 6.743 o-049 
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Table 2: Nucleus i^), bulge (^) and disk C") total magnitudes for HI Zw 2, 
Mrk 506 and Mrk 509. The errors of the magnitudes are shown as a superscript 
index to the corresponding values. The values of J^Wpsf and crgt are listed as 
well. 

Code /it^'ot /^L /^tot ^Wpsp (Tfit 

[ mag ] [ mag ] [ mag ] [ " ] [ / ^ " ] 



laS 


16.9350-113 


16.4350-193 


17.8310-028 


1.575 


0.025 


lay 


073 0.060 


15.7970-143 


16.2140-017 


1.560 


0.017 


lai?c 


^5 779 0.186 


15.2220-316 


15.5220-051 


1.311 


0.045 


la/c 


14.8930-121 


15.2270-341 


14.5120-023 


1.428 


0.047 


2bC/ 


15.8580-120 


15.4460-210 


I5.3530-053 


1.816 


0.013 


IbS 


16.591 




15.0460-010 


2.502 


0.028 


2bS 


16.725°-iii 


16.9560-244 


15.1580-023 


1.752 


0.035 


IhV 


153140.019 




14.2190-007 


2.350 


0.009 


2hV 


16.6040-082 


16.1500-107 


14.215 0-016 


1.516 


0.015 


IhRc 


16.230^-186 


16.2392-963 


13.6740-013 


2.360 


0.014 


2hRc 


15.8380-023 


15.8380043 


13.6410-010 


1.753 


0.007 


Ib/c 


15.6900037 


15.5540050 


12.9140-009 


2.293 


0.007 


2b7c 


15.8840037 


15.0770-034 


12.9570014 


1.459 


0.006 


3cU 


12.6640-074 


13.190 0-201 


15.4590-058 


2.711 


0.037 


IcB 


15.2880-063 


14.5410-086 


16. 134 0-045 


2.141 


0.020 


2cB 


15.1910-063 


14.2470-091 


I5.8I30-015 


1.861 


0.016 


3cB 


14.1400-145 


13.8690-212 


15.6960-022 


2.420 


0.024 


IcV 


14.5950035 


14.454 0066 


I5.OO4O-018 


1.877 


0.031 


2cV 


14.3980-034 


14.3810-082 


14. 799 0-027 


1.653 


0.022 


3cV 


14.169 0-076 


13.7930-144 


14.9010-028 


2.204 


0.014 


IcRc 


14.2020-020 


13.8740-034 


14.4320-011 


1.777 


0.017 


2cRc 


13.9200-022 


13.8230-049 


14.379 0-016 


1.491 


0.023 


3ci?c 


13.4240-023 


13.8160072 


14.2310-018 


2.252 


0.015 


Ida 


13.9290-020 


13.5600-030 


13.623 0-013 


1.797 


0.019 


2c7c 


13.6970-019 


13.7850-045 


13.5730-011 


1.977 


0.014 


3c/c 


13.2890-022 


13.5370057 


13.2160-018 


1.954 


0.017 
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Table 3: Bar parameters for III Zw 2. The bar lengths obtained from the ellipticity 
profiles employing the maximal and the minimal ellipticity criteria are listed in 
the last two columns for comparison with the bar length obtained from the SBPs 
decomposition. 



Code 


A'' cen 

[mag/D"] 


^ end 
["] 


m 


7 (max) 
^bar 

["] 


/ (min) 
'bar 

["] 


laB 


23.9410-043 


11.4300-142 


1.8740-110 


7.978 


11.350 


laV 


22.5100034 


11.2580-110 


1.85 2 0-091 


6.991 


11.214 


lai?c 


22.1890-111 


11.9060-™ 


2.5060.536 


6.726 


12.145 


la/c 


22.4090.140 


9.4680-213 


0.8760.190 


6.824 


11.252 



a sum of a Gaussian, a de Vaucouleurs and a truncated exponent was used as a 
model SEP. We have not found evidence of Freeman type II profile in Mrk 509 
to justify the usage of a truncated exponent. We have found a nearly exponential 
bulge, n ^ l,m all three galaxies decomposed (see also H^HHIl). 
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Table 4: Ring parameters for Mrk 506. 



Code 


cen 




^ cen 






[mag/D"] 




["] 


["] 


2hU 
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Ill Zw 2 - September 9/10, 1997 - V 
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Figure 1: SBP decomposition for III Zw 2. In the upper panel the observed 
SBP (open squares) and the decomposed profiles of the structural components are 
shown; the resulting model SBP closely follows the observed one. In the lower 
panel the residual profile is shown. 
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Mrk 506 - July 18/19, 1998 - 
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Figure 2: Tlie same as in Fig[T]but for Mrk 506. 



11 



Mrk 509 - July 20/21, 1998 - V 
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Figure 3: The same as in Fig[T]but for Mrk 509. 
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